New Insight into the Photoproduction of Z}*-meson in QCD 



A.V.Berezhnoy 

Nuclear Physics Institute of Moscow State University, 
Moscow, Russia 
and 

V.V.Kiselev, A.K.Likhoded 

State Research Center of Russia "Institute of High Energy Physics", 
Protvino, Moscow region, 142284, Russia 

Abstract 

In the framework of the model motivated by perturbative calculations in the 
fourth 0(afa)-order, the estimates for the 7p cross-section of the .D*-meson pro- 
duction in the ZEUS experiment are performed. We factorize the hadronization 
of (cg)-state hardly produced in perturbative QCD, which allows us to take into 
account the higher twist terms in the powers of 1/px at pt ~ ni c and to cor- 
rectly reproduce the c-quark fragmentation dominant at high pj< ^> m c to the 
given order of a s . We find a good agreement with the experimental data on the 
photoproduction, if the color-octet (cg)-state is taken into account, which yields 
(0 (8) } ss 0.33 - 0.49 GeV 3 . 

The measurement of the charm production by the ZEUS and HI collaborations at 
HERA in the interaction of quazireal photon with the proton |l|, Q rises the problem 
on the data interpretation in the framework of QCD. The description of process for the 
heavy quark production in the perturbative QCD (pQCD) is completely justified due 
to the hardness of partonic subprocesses. At present, there are two types of calculations 
based on the results in the next-to-leading order over a s , i.e. in 0{aa 2 s ): for the production 
of massive and massless quarks, correspondingly [[|, |J- 

In the both approaches the fragmentation model for the c-quark hadronization into 
D* was used at px 3> rn c . It was applied in the following form: 

d 2 a D * _ da C c{k T , jj) 
dzdpx dkx 

where D c ^e>*(z, fi) is a fragmentation function normalized by a probability of the c-quark 
fragmentation into the D*-meson, W(c — > D*), which is empirically determined from the 
e+e^-annihilation data 0, where W(c — > D*) = 0.22 ± 0.014 ± 0.14, and /x ~ p T is the 
factorization scale for the perturbative partonic cross section da/dhx- It is clear that 
such the approach does not allow one to describe the experimental data in the region 
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of pt ~ Tn c , where the fragmentation is not in the saturation regime, yet. The fact is 
valid even for the case of calculations with the massive quarks, wherein the pQCD can be 
correctly applied in the region of px ~ m c- 

The most important assumption of this method is that the factorization takes place 
for the probability of isolated c-quark hadronization and its hard production. This cir- 
cumstance automatically leads to the mentioned restriction: px 3> m c . 

In this paper, we discuss a model of the charmed meson production with the hadroniza- 
tion of (cg)-system into the D*-meson instead of the hadronization of isolated c-quark. 
So, it is quite reasonable to apply this approach in the region of px ~ m c , too, wherein the 
pQCD can be used. In the production of (eg) -state, it is necessary to take into account 
the production of the additional quark pair, so that, in the leading order of a s , the 0{aa z s ) 
calculation must be used. Among the full set of the tree level diagrams on Fig. 1, we 
can see several diagrams (i.e. 16 and 19) of heavy quark fragmentation, which dominates 
at pt 3> m c . At high transverse momenta, they reproduce the results of fragmentation 
model with the following fragmentation function ||: 

D *-°- w = ^F(iXo> - 2(3 - 2r)z + 3(3 - 2r + 4rV 

-2(1 - r)(4 - r + 2r 2 )z 3 + (1 - r) 2 (3 - 2r + 2r 2 )/), (2) 
where r = m q j (m q + m c ) and 

(O eff ) « (0(i)> = ^ + ^f) (D*(p)\(^q)(q ll/C )\D*(p)}, 

which corresponds to the squared value of the quarkonium wave function at the origin for 
two quarks in the framework of nonrelativistic potential model: (0(i))|at.r = 1^(0) | 2 . We 
have introduced the effective mass of light quark m q , which determines both the normal- 
ization and the form of fragmentation function. It is worth to stress that, in the higher 
orders of QCD, the hard part of processes does not depend on m q . However, the problem 
of infrared divergences always presents in these cases. Usually, it can be removed by a 
redefinition of initial and final states. For the initial states one takes into account the 
Q CD-evolution of intrinsic quark-gluon sea. One can introduce a regularization parameter 
(/it in the NLO-calculations in (§) to avoid the difficulties connected to the nonper- 
turbative hadronic sea in the final state. In the model under consideration, m q ~ m p /2 
is used as an infrared regulator. In this case the gluonic virtualities are greater than 
\kg\ ~ p T + m 2 p . So, the fact allows us to believe in the reliability of pQCD calculation. 
The Z}*-meson production in e + e~-annihilation has been used to fix the model parame- 
ters. We find that the form and normalization of fragmentation function with the account 
for the QCD-evolution is in rather a good agreement with the LEP data, if the following 
values of parameters are chosen as below: 





= 0.35, 


m q 


= 0.3 GeV, 


m c 


= 1.5 GeV, 


ie//\ 


= 0.25 GeV 3 


D*) 


= 0.22. 
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In e + e -annihilation the enforcement of color-singlet (cg)-state takes place in compar- 
ison with the octet contribution due to the color amplitude structure, so that 

(O eff } = (0 {1) }+ 1 -(0 {8) } } (4) 

where 

(0( 8) ) = M + ^) (D*(p)\(n»\ a q)(qiA b c)\D*(p)) 

It is qualitatively clear that the relative contribution of the octet four quark operator 
(0(g))/ is determined by the relative velocity v of the quark motion in the hadron. 
For example, in the hard production of heavy quarkonium the octet fraction is strongly 
suppressed due to the nonrelativistic expansion in v <C 1 0. For the meson, containing 
the light relativistic quark, we can expect that (0( 8 ))/ (0(i)) ~ 1, since v ~ 1. Supposing 
this arrangement in e + e~-annihilation, we see that (0(i)) is dominant in (O e ^). The 
normalization of octet operator (0(8)) remains the free parameter. It can be determined 
from a process with another color structure of the production amplitude, for example, 
from the photoproduction. In this case, the essential contribution due to diagrams of 
nonfragmentation type (recombination) takes place at ~ m c . We have to emphasize 
that, in this kinematical region, the octet contribution is not suppressed in comparison 
with the singlet one. 

The calculation results are shown in Fig. 2 in comparison with the experimental data 
of ZEUS in the following kinematical region: pt > 2 GeV, — 1.5 < 77 < 1.5, 130 GeV < 
W < 280 GeV, Q 2 < 1 GeV 2 , where W is the total energy of 72? interaction, Q 2 is the 
photonic virtuality, and rj is the pseudorapidity of Z}*-meson. 

We see that the value of the D*-meson production cross section cannot be explained by 
the singlet mechanism only. The relative value of octet contribution (0(8))/(0(i)) = 1.3 
allows us to describe the spectra of D*-mesons. We would like to emphasize that the 
singlet term, say, at px > 4 GeV is very close to the predictions in the model for the 
production of isolated c-quark |3|, [|] . 

Following a request by the ZEUS Collaboration we have calculated the differential cross 
sections for the Z)*-meson production, as predicted by the model under consideration and 
shown in Fig. 3, at low W for the following kinematical region: pt > 2 GeV, —1. < i] < 1.5, 
80 GeV < W < 120GeV, Q 2 < 0.01 GeV 2 . The same values of (0 (1) ) and (0 (8) ) have 
been used in the both calculations. 

Using the experimental data by the ZEUS collaboration, only, we cannot exclude the 
possibility of other ratio between the singlet and octet operators, if the octet fraction is 
not suppressed in the fragmentation regime in e + e~-annihilation. Then, for example, at 
(O(i)) ~ 0.14 GeV 3 and (0( 8 ))/ (O(i)) ~ 3.5 we find a good description of the experimental 
data (see Fig. 4). However, the low value of singlet operator in combination with the high 
fraction of octet contribution results in quite a slow variation of the octet term. 

Further, in the calculations we have neglected the differences between the spectra 
of color-octet vector and pseudoscalar states, whose contributions have been effectively 
summed up in the corresponding operator. Moreover, we have not considered a possible 
influence by the octet P-wave states, which seem to be very unstable for the heavy-light 
system. 
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In conclusion, we have presented the model for the charmed meson production in 
photon-proton collisions with the hadronization of (cg)-state in contrast to the usual 
restriction by the isolated c-quark. It allows us to make justified predictions in the range 
of pt ~ m c , which improves the results of simple c-quark fragmentation regime dominant 
at high transverse momenta. In addition, the model includes the contributions of color- 
singlet and octet (cg)-states, which depend on the color structure of production amplitude. 
So, in e + e~-annihilation the octet fraction is suppressed, while in the photoproduction we 
have found that the singlet contribution underestimates the data, and the octet operator 
can be measured so that, in the conservative evaluation, we extract 

(0 (8) ) « 0.33 - 0.49 GeV 3 . 

To improve the estimation, we plan to analyze the model predictions for the electro- 
production, where the photon can be essentially virtual, in comparison with the current 
data. 
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Fig. 1. The leading order diagrams for the (cg)-state 
production in #7 interactions. 
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Fig. 2. The distributions of D* ± -meson photoproduction for 
(0$)/{Oi) = 1.3 in comparison with the ZEUS data 
at 130 GeV < W < 280 GeV and Q 2 < 1 GeV 2 . The 
color-singlet term is shown by the dashed histogram, the 
octet contribution is given by the dotted histogram, and 
the solid histogram presents the sum of two terms. 
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Fig. 3. The distributions of D*-meson photoproduction for 
(0 8 )/(Oi) = 1.3 at 80 GeV < W < 120 GeV and 
Q 2 < 0.01 GeV 2 . The notations are the same as in 
Fig. 2. 
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Fig. 4. The distributions of D*-meson photoproduction for 
(0$)/{Oi) = 3.5 in comparison with the ZEUS data 
at 130 GeV < W < 280 GeV and Q 2 < 1 GeV 2 . The 
notations are the same as in Fig. 2. 
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Fig. 5. The distributions of D*-meson photoproduction for 
(0 8 )/(Oi) = 3.5 at 80 GeV < W < 120 GeV and 
Q 2 < 0.01 GeV 2 . The notations are the same as in 
Fig. 2. 
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